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Abstract 



A systematic search for different viable models of the dark energy universe, all of which 
give rise to finite-time, future singularities, is undertaken, with the purpose to try to find 
, a solution to this common problem. After some work, a universal procedure to cure all 

future singularities is developed and carefully tested with the help of explicit examples cor- 
responding to each one of the four different types of possible singularities, as classified in the 
literature. The cases of a fluid with an equation of state which depends on some parameter, 
^ |. of modified gravity non-minimally coupled to a matter Lagrangian, of non-local gravity, and 

D ' of isotropic turbulence in a dark fluid universe theory are investigated in detail. 

^ " 

PACS numbers: 98.80.-k,04.50.-f h 

> 

CN : 1 Introduction 

One of the most important problems of modern cosmology is to explain in a plausible way the late-time 
, acceleration of the universe expansion. Recent observational data, coming from more and more precise 

' high-redshift surveys of type la supernovae P, Q and from anisotropics in the power spectrum of the 

' I , cosmic microwave background 0, |j, indicate that, right now, our Universe is undergoing a phase of 

accelerated expansion. The favored explanation for this behavior is that the universe is filled with some 
form of dark energy (DE). However, none of the DE models proposed so far is completely satisfactory. 
In order to explain the late-time cosmic acceleration, and to unify it with inflation, a new family of 
theories has been proposed, based on modifications of Einstein's gravity (see Q for a recent review). 
These theories, such as f{R) and /(G) modified gravities (R for the curvature and G for the Gauss- 
Bonnet invariant), and also non-local gravity, have the power to unify primordial inflation with late-time 
inflation in a natural way. While a number of the original models had to be rejected by several reasons, 
some subfamilies of them (see 0) fulfill all available cosmological constraints and have been checked to 
pass the solar system tests and cosmological bounds, e.g., the fact that, at these scales, Einstein's gravity 
is valid to a very high degree of accuracy. But even these, so called viable, modified gravities are not 
free from other problems, one of the most important being the frequent appearance of finite-time, future 
singularities. 

The singularity problem is indeed of fundamental importance in modern cosmology. In order to 
address this issue rigorously it would be necessary to develop a fully-hedged theory of quantum gravity, 
but this has proven to be a very difficult, up to now impossible, challenge. Anyhow, the presence of a 
finite-time, future singularity may cause various problems of physical nature, as instabilities in current 
black hole and stellar astrophysics. It turns out that, even without the support of a quantum gravity 
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theory, it is still meaningful at first instance (and of great importance) to try to find natural scenarios, 
already at the classical and semiclassical levels, that may cure the emergence of these finite-time, future 
singularities. Usually, one starts with some given theory and then solves the corresponding equations 
of motion in order to define the associated background dynamics. But, for several models of modified 
gravity or scalar-tensor gravity, there is in fact another possibility owing to the fact that those models 
are defined in terms of some arbitrary functions or potentials. The new possibility consists in using 
the freedom in the choice of such arbitrary functions or potentials, with the aim to reconstruct the — in 
general complicated — background cosmology which complies with the latest observational data. 

In this paper, a reconstruction program of this sort is applied to a number of theories which do give 
rise to models that exhibit finite-time, future singularities, with the purpose to investigate its structure 
in detail and try then to cure this common problem. Specifically, it will be found in all cases that the 
addition of an term provides a universal tool capable to cure these finite-time future singularities. 

The work is organized as follows. In Sect. 2 the case of a fluid with an equation of state (EoS) 
depending on a parameter, a, which can give rise to finite-time future singularities, is considered. Different 
possibilities for the evolution of this universe, depending on the value of a, were studied in [§]. In the 
subsections it will be explicitly shown that, adding a function G{H, H, H, ...) to the EoS of the fiuid, the 
different singularities can indeed be cured. Such function G is actually to be considered as a modification 
of Einsteinian gravity. Sect. 3 is devoted to the case of non-minimal coupling of modified gravity to 
a matter Lagrangian. The reconstruction scheme is developed for this case and the specific example 
of the cosmology given by the Hubble function H{t) = hs/{ts — t) is analyzed. The calculation of the 
Friedmann equations for this example of modified gravity non-minimally coupled to the matter Lagrangian 
is explicitly carried out in App. A. In Sect. 4 the case of non-local gravity is discussed. The example 
of the de Sitter space is reproduced in the framework of non-local gravity. It is pointed out there that 
theories of this kind can give rise to a finite-time, future singularity. Finally, Sect. 5 is devoted to the 
case of isotropic turbulence in the dark fiuid universe. It will be shown there that the contribution of 
the turbulent part of dark energy can be reproduced through the use of a scalar-tensor theory. Several 
examples are discussed in detail. The paper ends with some conclusions and an outlook. 

2 Accelerating universe with and without a future singularity 

Different accelerating universes, with and without finite-time future singularities, are considered in this 
section. We work with a particular fiuid with EoS given by: 

P = -P + Ap", (1) 

The nature of each singularity depends on the value of the parameter a. All possibilities, corresponding 
to the different values of a, have been studied in [8]. It will be shown below that introducing a specific 
function, G{H), into Eq. ((T|), 

p^-p + Ap'- + G{H), (2) 

the singularity can be avoided. 

The spatially fiat Friedmann-Robertson- Walker (FRW) universe in the frame of General Relativity 
will be considered 

ds^ = -dt^ + a{tf (dx^ + dy^ + dz^) , (3) 
a{t) being the scale factor. The Friedmann equations are 

= yp, (4) 

i^ + i?' = -y(p + 3p), (5) 

with = SttG. For future use, it is useful to classify the future singularities as in namely, 

• Type I ("Big Rip") : For t — t^, a — oo, p — cx) and \p\ — > oo. This type of singularity is discussed 
in 0. 

• Type II ("sudden") : For t ^ tg, a ^ as, p ^ Ps and \p\ oo. 
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• Type III : For t t^, a ^ as, s-oo and \p\ — ^ oo. 

• Type IV : For t ^ ts, a ^ Us, p ^ Ps, \p\ ^ Ps and higher derivatives of H diverge. 

Here tg, as, Ps and Ps are constants, with as ^ 0. 

In the next subsections we wih consider four different cases for the dark fluid H]), which lead to the 
four possible different singularities. We will produce a particular function G{H) for each case, which will 
cure each specific singularity. Finally, in the last subsection, a general function G{H, H, H, ...) that cures 
all possible finite-time future singularities of (HJ will be constructed. 



2.1 p=-p+Ap^ 

For this particular EoS a Type III singularity occurs {H{t) oc {to — t)^-^" , a > 1). However, if the 
following EoS is considered 

p=-p + Ap^ + G{H), (6) 

with 

GiH)^-^H^ + C, (7) 

being C real, then the singularity is avoided. In order to explain this fact, one must take into account 
that, using the Friedmann equation ([4]), Eq. (O reduces to 

GiH) = -Ap' + G, (8) 

and using now Eq. ([5]) the EoS ([5]) yields 

P^-P + C, (9) 

which does not have future finite-time singularities of any kind. It is interesting to note that the above 
specific choice of G{H) can be motivated by modified gravity (see Q). 

2.2 p=-p + Apl 

_ 2 

In this case, a Type I singularity occurs, namely H{t) — j^^, a = 1 and A < Q or H{t) oc 
(to ~ t)'^^ , 1/2 < a < 1. Considering now Eq. with 

G{H)=BH^, (10) 

and using the Friedmann equation given by the EoS © reduces to 

p ^ -p + Ap^ + B'p, (11) 

where B' = B^. When p ^ oo the term B'p dominates over the term Api which was the one that 
caused the singularity. If B > then this Type I singularity is removed. 

2.3 p=-p+Ap^ 

With this EoS there appears a Type IV singularity, namely H{t) cx {to — t) ^^^^ and i_2a ^"^^ integer. 
Now, if Eq. © is considered, with 

G{H) = C, (12) 

where G is real, then the EoS ^ turns into 

p = -p + Ap-^ +G. (13) 

Whenever p — !■ the term G dominates over the term Aps and the existing singularity, of Type IV, is 
cured. 
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2.4 p = -p + Ap~'^ 

In this case a Type II singularity shows up, namely H{t) oc (to — t) , a < 0. By taking in Eq. ^ 

G{H) = -A^H-^ + C, (14) 

with C real, and using Friedmann's equation Eq. (IT4l) yields 

G{H) = + c. (15) 

Then, the EoS (O reduces to 

p=-p+C, (16) 

which does no more exhibit any kind of finite-time future singularity. 

Summing up, we have shown in all previous situations that, for each case, a particular function, 
G{H), can be found which cures the singularity which can possibly appear in the model given by Eq. ([1}. 
It is interesting to realize that a function G{H) of this kind can be interpreted as a contribution of 
modified gravity, as was shown in [^, [llj. Putting everything together, we thus have demonstrated, in a 
very explicit way, how modified gravity is able to cure all finite-time future singularities that can possibly 
appear in a fiuid with the particular EoS given by Eq. ([l}. 



2.5 A generic function G{H, H, H, ...) which cures all singularities for the fluid 
with p = —p + Ap"' 

In this subsection we will develop a systematic method for finding a function G{H, H, H, ...) which avoids 
any possible singularity for the model with p = — p + Ap°'. Let us recall [ll[ that every F(i?)-modified 
gravity can be seen as Einsteinian gravity with a particular EoS which absorbs the effects of F{R). 
We consider the case 

F(i?) = i? + /(i?), (17) 

being f{R) = aB? (it is known that the term E? cures the possible appearance of all future finite-time 
singularities [1, [l2, 13|). For Eq. P7)) . using the results obtained in fll], it follows that 



-\f{R) + 3 + i/j f'{R) - 18 {ah'^H + HH^ f"{R) 



Pmatter ? 



i/(i?) - (SH^ + f'{R) + 6 {d,H'^H + 4:H^ + 6HH + /"(i?) + 



+ 36[AHH + Hj f"'(R) +PraaUer- 

It is also known that the Friedmann equations can be written as 



(18) 



(19) 



(20) 



P^ff = 



2H + 'iH' 



(21) 



In the case of f{R) — aR^ and taking into account Eqs. (|20l) and (l?ll) . Eqs. ((T5)) and (|2.5p reduce to 

18a 



18a / • 

Pmatter — Peff J" I ^ 

6a 



- 6H^H - 2HH 



Pmatter = Pef f - ^ {^H^ + ISH^H + l^HH + 27j) , 

respectively. If we now consider the EoS 

Pmatter — Pmatter ~^ -^Pmatter ^ 



(22) 
(23) 

(24) 
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introducing into Eq. (j24L the results obtained in Eqs. ((22)) and (j23L we get 

Peff = -Peff + Ap-^ff + G{H, H, ...), 

where G{H, H, ...), in the case of F{R) = R + aR^, is given by: 



(25) 



A 



18a ( - GH'^H - 2HH 



"-(3i72)"|. (26) 



Thus, using F{R) modified gravity (adding in the action a term proportional to R^, see [1, [l2, 13|), a 
function G{H,H, ...) has been found which cures all the singularities which appeared in the model given 
by Eq. (j24p . in the frame of Einstein's gravity. 



3 f{R) Modified gravity with possible future singularities for 
the case: jC = + f{R)Cm 

In this section we will investigate f{R) modified gravities non-minimally coupled to matter-like La- 
grangians that lead to future finite-time singularities. A general Lagrangian density of this sort is (see 

C=^R+fiR)Crn. (27) 

By varying with respect to the metric, the following field equations are obtained 

^ (r^^ - + R^uf'imCr,, + (.9p.V2 - V^V,) {f'{R)C,n) - lf{R)Tf,u - 0. (28) 

Furthermore, considering the particular case given by C„i = d^4'd^4'^ and varying with respect to the 
field, we get 

{^f{R)d^^<f) = 0, (29) 
and, if it is now assumed that (j) — (l>{t), Eq. reduces to 

dt {^9f{R)d'^) = ^ ^f{R)M - 0. (30) 

Considering a spatially fiat FRW universe, this is = a{t)^ , and writing dt4> — 4>, Eq. ([50)1 yields 

a{tff{R)^^C, (31) 

where C is a constant. Then, 

Taking once more into account that cj) = 0(i), one gets 

^" " " 'a{tff{Ry ^^^^ 

Thus, considering a spatially- flat FRW universe and Eq. Friedmann's equations for the Lagrangian 
density given by (l?fl) can be derived. Details of the long calculations leading to these equations are given 
in appendix [XI One should note that (|X]) and constitute a pair of differential equations for f{R) 
with the variable being the scalar curvature, R. Thus, starting from a given Hubble function H{t) and 
taking into account that t = t{R), from the relation R — 6H{t) + 12H{t)^ and by using (jX]) or (|X]) we 
obtain a specific function f{R) that reproduces the given Hubble function. 

Another important remark is the following. From the stress-energy tensor T^^, 
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by comparison of 



with ((M)) . we obtain the result: 



2 5{^C,, 



Too = -p.goo = P 
Til = pgii = a{t)'^p 



a{tYf{Rf 



(35) 



(36) 



Now, from ([551) we know that, for this model, a singularity of type II is avoided (for the classification of 
the future singularities, see Q or Sect. II above). 

In order to find a model with the Lagrangian density (|27D that ends in a future finite-time singularity, 
we now consider the case: 

H{t) = 7^ 



where, from R = 6H{t) + 12H(tf, it follows that 



i = <.s - 



hi 



R 



where hi = ^6/is(l + 2hs 



Then, we can write 
H 



h2R^, 

H ^ hsR^, 
H = hiR^, 



(37) 
(38) 

(39) 

(40) 
(41) 



being h2 



&{l+2hs 



and /i4 = 



t{R) 



H{t')dt' = In 



3{l+2hs)^y6h,{l+2h,) 
ts — to 



For ((37| . we have 



hi 



(42) 



to 



hence, Eq. 



J reduces to 
f{R) + aiR'^-f{Rf 



asR 



df{R) 
dR 



2a4 



i?2 



fdf{R)\ 
f{R) \ dR J 



(43) 



where ai = -^§^4 j , 02 = 3/is + 1, 03 = /13 + T/ij and 04 = 6/i2(/i4 + 4/12/13)- The solution for 

Eq. gives us the function /(i?) that reproduces the Hubble function given by (1571) . Eq. (H5)) can be 
solved, but the solution obtained is too long and not particularly insightful to be written here. 

To conclude, we should mention that quantum gravity effects (which usually contain different powers 
of the curvature) become very important near the future singularity (see [IBl)- There, classical consider- 
ations are not valid. It is known [l^, [l^] that the DR term works against the singularity. Thus, an R^ 
term (which will generate a OR term) would cure the possible singularities that could arise in the theory. 



4 Future finite-time singularities in non-local gravity 

The case of non-local gravity will be here considered. This theory gives a natural unification of inflation 
with the current cosmic acceleration and it is inspired by quantum loop corrections (see 0, 17, isl. [l9j). 
Non-local effects come from the introduction in the action of the inverse of the D'Alembertian, 0~^, and 
the simplest action of non-local gravity is therefore 



S = Jd^xV^ 1^ [R{l + f{D-^R))-2A] +£matter{Q;g)y 



(44) 
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where Q stands for the matter fields and A is the cosmological constant. Introducing two scalar fields, ry 
and ^, action can be rewritten as 

^d^xV^ |^[i?(l + /(ry))+e(D??-i?)-2A]+£™,„er| (45) 

If we assume a spatially-flat FRW metric, and that ry = ri{t) and ^ = ^(t), the equations of motion for 
the scalar fields and the Friedmann equations read (see Q) 

Q = i) + ZHi] + QH + 12H^, (46) 

Q = i + mi-m+l2H^)f[7^), (47) 

= -iH^ (1 + /(t?) -0 + \iv - (/'(7y)7) - e) + A + n^Pmatter, (48) 

0=(2H + 3H^) (1 + fiv) - + ^^^7 + + 2^^^) (/(^) - - A + «;'p™atte.. (49) 

Given a Hubble function, H{t), Eq. (|46|) can then be solved to obtain rj — ri{t); moreover, the function 
^ = £,{t) can be obtained from Eq. (|Tf)) if one assumes a form for the function /{rj). Once we have the 
functions r]{t) and ^(t), assuming an EoS for Pmatter and Pmatter, Eqs. P5|) and (|49l) yield the relation 
between the different parameters that appear in the model (i.e. the constants of integration of the 
functions ri{t) and £,{t), the cosmological constant A, etc.). 

In the previous scheme is used to show that de Sitter space {H{t) — Hq) can be a solution in 
non-local gravity. This is the case for matter of constant EoS w, when 

H^Ho, (50) 

Tjit) = -4Hot, (51) 

fiv) = foei=foe-^, (52) 

with Hq, /o, /3 and being constants. For de Sitter space and for matter with constant EoS uj, the 
energy density is 

Pmatterit)^P0e-^<^'+^^"°'. (54) 

In order to fulfih Eqs. (03, (HZl), (011) and (gS]), it is necessary that 

^" ^ si/foTsZ;) ' (^^^ 

6=^-1 + ^. (57) 

Thus, de Sitter space can indeed be a solution of non-local gravity. 

There are also singular solutions; however, they are very involved and will not be considered here. 
These singular solutions could also be cured by the addition of an term (inspired by quantum gravity 
effects near the singularity), a procedure which could a gain turn into a universal tool in order to suppress 
all finite-time future singularities, as before (see (gI [i^Ti3|). 
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5 Reproducing isotropic turbulence in a dark fluid universe with 
scalar-tensor gravity 

We will now emphasize the fact that a scalar-tensor theory can be used \20\ in order to reproduce isotropic 
turbulence in a dark fluid universe 2l|. To this end, let us consider the following scalar-tensor theory 
action 



which leads to the Friedmann equations 



(58) 



= -f V{(i>) + Pmatter, (59) 

-y{<i>) ^ V^atter. (60) 

If we consider that the scalar part of the action dominates over the matter part, these Friedmann equations 
reduce to 

= lc.{^)^^ + F(0), (61) 

2,H -\- 311^ 1 , -o , ^ , , 

-w(0)02 _ v{<j,). (62) 

On the other hand, the Friedmann equations describing isotropic turbulence in a dark fluid universe (see 
2l|) are 

^ ~ Pdark H~ Pturb ~l~ Prad ~^ P matter i (63) 

2H + 3i/2 

5 = Pdark + Pturb + Prad + Pmatter ■ (64) 

Now, under the proviso that the turbulent part pturb dominates, Eqs. (j63p and (j64p reduce to, respectively. 



2H + 3i?2 



Pturb, (65) 



Pturh- (66) 



Thus, comparing Eqs. §1^ and (162]) with Eqs. (|65j) and ([66 

1 



2 

and assuming at this point that cf) — t, we get the expressions 



Pturb = §^(<^)0' + Vicf>), (67) 
= luji^H^ ~ V^(</)), (68) 



Pturh = ^C^(0) + V^(0), (69) 

Pturb = ^C^(0) - "^^(0), (70) 

respectively. Finally, the functions uj{(f>) and V{(p) are given by 

^{4') = Pturb + Pturb, (71) 

^(0) = Pturb - Pturb- (72) 

In '^l'] some examples of isotropic turbulence were given. In the following subsections, by using 
Eqs. dTlT) and (|72L scalar-tensor gravities that reproduce these specific examples will be constructed. In 
what follows below, the EoS pturb = ijJturbPturb is assumed for the isotropic turbulence. 



8 



5.1 De Sitter space 

In this subsection we consider the case of de Sitter space, for which the scale factor is given by a{t) 
aoe^°*, with constant oq and Hq. For this case, the energy density of the turbulent part is 



Pturb — G 



c 



turb 



2Cna, 



0«o 



57; 



turb 



(73) 



with C — [pturbitin)] ^ and where tin is an initial time, from which the universe starts its develop- 
ment onwards. Cq is an integration constant, which can be determined from the initial condition t = tin 
and jturb = 1 + ^jJturb- Then, use of Eqs. ([7T|) and ([7^ yields 



Jturb e 



-Sjtui-bHoi 



^{4>) = Pturb + Pturb — IturbPturb 

c 



e 2 



IturbHati. 



e 2 



2Coa| 



T^turb^^O 







and 



(1 - ioturb) e-3^--^«* 



Sjturb 

^(</') — Pturb — Pturb = (1 ^ ^turb)pturb 

c 



^Iturb 



(74) 



37> 



' g-^lt^LTbHoti-a _ ^-%-ft^rbHo', 



turb 



IturbHo 



57*1 



urb 



(75) 



As a consequence, considering the action ([55]) . with Eqs. (j5.ip and (15. ip . together with the assumption 
of being (f) = t, is just equivalent to consider isotropic turbulence in a de Sitter universe. 

5.2 Effective quintessence-like power-law expansion 

In this subsection we discuss the situation when the scale factor is given by a{t) — aQt^° , with gq and Hq 
constant. For this case, the energy density of the turbulent part is 



pturb — <- 



5C 



6 (l - ^^turbh, 



o) ^ ^ 



(76) 



being C = [pturb{tin)] ° and tin, as before, an initial time from which the universe starts its 
development onwards. Again, Co is a constant of integration, to be determined from the initial condition 
t = tin and 7tur& = 1 +<jJturb- Repeating the procedure of the previous subsection, we can write 



Jturb ^ 



^{4>) = Pturb + Pturb — IturbPturb 

5C 



and 



(1 - LOturb) i-37*„..'.o 



(A-h*^.bh0 _^-hturbha\ „ f 

6 (1 - ^Iturbho) V / 
^(</') = Pturb — Pturb = (1 ^ ^turb)pturb = 

^'^ (fl-ht^.bho _ A-hturbho\ ^ f 

6 (1 - ^Jturbha) V / 



Iturbho 



(77) 



Iturbho 



(78) 



Also, as remarked above, the scalar-tensor gravity given by the action ([55]) . with Eqs. (|5.2I) and (|5.2p . 
together with the assumption cf) = t,\s here equivalent to consider isotropic turbulence in the dark energy 
component, now in an effective quintessence-like power-like expanding universe. 
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5.3 Phantom-like power law expansion 

Here we will deal with a phantom- like power law expansion, given by a{t) = ao{ts — , which is known 
to give rise to Type I finite-time future singularities, with constant oq and ho. The energy density of the 
turbulent part is 



5C 



6(1 + ^Jturbho) 



Iturbho 



(79) 



where C = [pturb{iin)\ and tin and Co have exactly the same meaning as before. In this case 



'^{4>) = Pturb + Pturb — IturbPturb 

5C 



6 (l + f Ttnrb/lo) 



{ts - ^--''O - {ts - t,„)l+2^--''« 



(80) 



and 



{1 — LOturb) {ts~t) 



3-(turbho 



^{'l>) — Pturb — Pturb — (1 ^ '^turb)pturb 

5C 



6 (l + \lturbho) 



{ts - t) 



{ts ^m) 



0% 



\ltiirbho 



(81) 

Therefore, a scalar-tensor gravity with functions uj{<j}) and 1^(0) given by expressions (j5.3p and (j5.3p . 
respectively, is equivalent to isotropic turbulence in a phantom-like power-like expanding universe. 

Thus, we have shown that an equivalence exists between isotropic turbulence and scalar-tensor gravity 
and that in the isotropic turbulence theory, finite-time future singularities appear too (see Sect. 15. 3p . As 
in previous sections, addition of a term to the Langrangian density can avoid the development of these 
finite-time future singularities [12, [l3| ■ 



6 Conclusions 

In this paper, a reconstruction program has been dealt with which uses the freedom in the choice of 
arbitrary functions or potentials, for several models of modified gravity or scalar-tensor gravity, with 
the aim to reconstruct a background cosmology — quite complicated in general — which complies with the 
latest observational data. Along this line, a systematic search for different viable models of the dark 
energy universe, all of which give rise to finite-time, future singularities, has been undertaken, having 
as goal their detailed study to try to find common features, in the search for a general solution to this 
important problem. Specifically, it has been checked that the addition of an term provides indeed a 
universal tool to cure these finite-time future singularities. 

More specifically, a universal procedure to cure all future singularities has been defined and carefully 
tested with the help of explicit examples, corresponding to each of the four different types of possible 
singularities, as classified in the literature. To start, the case of a fluid with an EoS which depends on a 
parameter a, which can give rise to finite-time future singularities, has been considered. We have shown 
explicitly that, adding a specific function G{H, H, H, ...) to the EoS of the fluid, the different singularities 
can be cured, and it has been seen that this function can actually be considered as a modification of 
Einsteinian gravity. 

The case of the non-minimal coupling of modified gravity to a matter Lagrangian has been investi- 
gated, too. The reconstruction scheme has been run for this case and the example of the cosmology given 
by the Hubble function H{t) — hs/{ts — t) was analyzed. In Appendix A, the calculation has been done 
of the Friedmann equations for this non-minimal coupling of modified gravity to a matter Lagrangian. 
Further, the case of non-local gravity has been discussed. The example of the de Sitter space has been 
reproduced in the framework of non-local gravity, having been pointed out that such kind of theories can 
also give rise to finite-time, future singularities. Finally, the case of isotropic turbulence in the dark fluid 
universe was discussed, with the conclusion that the contribution of the turbulent part of dark energy can 
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indeed be reproduced through the use of a scalar-tensor theory. As aheady indicated, several examples, 
corresponding to the different cases, have been studied in the paper in detail. 

Concerning future perspectives, it is rather clear that, in order to address the singularity issue in 
all rigor it will be necessary to develop a fully-fledged theory of quantum gravity, what has proven 
to be a very difficult, up to now impossible, task. In any case, the presence of a finite-time, future 
singularity may cause various problems of physical nature, as instabilities in current black hole and 
stellar astrophysics. And, even without the recourse to a quantum theory of gravity, it is still meaningful 
to try to find natural scenarios, already at the classical level, that may cure this possible finite-time, 
future singularities. This has been successfully addressed in the paper, with the explicit construction of 
a general, universal procedure to cure all future singularities — what has been carefully tested with the 
help of specific examples corresponding to each one of the four different types of possible singularities, as 
classified in the literature. 

To conclude, we should mention that quantum gravity effects (which usually contain different powers 
of the curvature) may become very important near future singularities. Even if classical considerations 
are, in principle, not valid there, it is known that the OR term works against the singularity. Thus, an 

term (which will on its turn generate a \3R term) would, in principle, cure the possible singularities 
that could arise in the quantum theory, too. As a consequence, our method here could presumably be 
also extended to the quantum case without much trouble. 
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A Friedmann equations for the Lagrangian density (1271) 



In this appendix we carry out, in some detail, the calculations which are necessary to obtain the Friedmann 
equations for the Lagrangian density of Eq. (P7|) . Assuming a spatially-flat FRW universe and taking 
into account Eqs. ^ and Eq. {t,t) from Eq. ^ reads 



3i7(i)2 



(ij(i) + H{tf 



aitr/iRy 



f{R) 



a{tff{Ry 



1 C2 



and Eq. (i,«) from Eq. ((28)) . 

2H{t) + iH{tf 



C {H{t)+3H{ty 



'Cdtdt 



f{R) 



1 C2 



- 2C^H{t)dt 



2a{tYf{R) 

f'jR) 
ait)<^f{Ry 



= 0, 



AtrfiRr 

Taking derivatives with respect to time, one gets 



2a(t)6/(i?) 



= 0. 



f'iR) 



and 



dtd, 



a{tff{RY 
fiR) 



1 



a{tff{Ry 



riR)R-6f{R)H{t)^ 



2f'{RfR 

fiR) 



a{tff{Rf 
+f"{R)R? + 2!\Rf 



a{tff{RY 

( l2RH{t) - R 



6/'(i?) {QH{tf - H{t)] + /"(i?) [R - \2RH{t) 



f'iRf 



i?2 



f{Rf 



f'{R)f"{R) 



fiR) 



Introducing Eq. ([M)) into Eq. (15^ and taking into account that 



a{t) = ao exp 



H{t')dt' 



(82) 



(83) 



(84) 



(85) 
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yields 



alH{ty exp 



:fiR) 



6 / H{t')dt' 

to 



f{RY + [H{t) + 7H{ty]f'{R) + 



^-^fiRf RH{t)r{R) = 0, 



(86) 



and putting Eqs. ((Ml) and © into Eq. this reduces to 



2^^^) + ^2^2 «oexp 



6 / H{t')dt' 



f{RY + 7{Hit)-3H{tr\nR) 



R 10RHit)\ + ( lORHit) - R] /"(i?) + 

i?2 



7iR) J 



6-^/'(i?)/"(i?)-i?Y"(i?) = 0. 



(87) 



The other possibihties for Eq. ([28]) are identities. 

We also know that R — 6H{t) + 12H{tY, which could be solved in terms of i as i = t{R). Taking 
this into account, Eq. (f5| can be written as 



af,H{t{R)fexp 



:f{R) 



t{R) 

6 / H{t')dt' 

to 



■f{Rf + H{t{R)) + 7H{t{R)f 



df{R) 
dR 



l2H{t{R)) H{t{R)) + AH{t{R))H{t{R)) 



+ ■ 



f{R) 

H{t{R))+AH{t{R))H{t{R)) \ H{t{R)) 



(df{R)\ 
\ dR J 

dy{R) 

dR^ 



0, 



and Eq. ^ as 



1 2H{t{R))+3H{t{R)r 
p(R) + aoexp 



tiR) 
6 J H{t')dt' 

to 



!{Rf + 7 H{t{R)) - 3H{t{R)f 



dfjR) 
dR 



^ -AOH{t{R)yH{t{R)) + m{t{R)f - 6H{t{R))H{t{R)) + HjtjR)) f dfjRW 

f{R) [ dR ) 



216- 



HitiR))+AH{t{R))H{t{R)) 



f{Ry 



[df{R)Y 
\ dR J 



-6 



-mH{t{R)fH{t{R)) + m{t{R)f ~ 6H{t{R))H{t{R)) + Hit{R)) 



d^fjR) 
dR^ 



216- 



H{t{R)) + 4H{t{R))H{t{R)) 



df{R) d^f{R) 



fiR) 



dR dR^ 



36 



(i(i?)) + AH{t{R))H{t{R)) 



d'fjR) 



= 0. 

(89) 

Eqs. (j^ and (j^ are the Friedmann equations for the Lagrangian density given by ^7} and constitute 
the two differential equations we were looking for f{R). 
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